Strong-correlation effects in Born effective charges. 
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Large values of Born effective charges are generally considered as reliable indicators of the genuine 
tendency of an insulator towards ferroelectric instability. However, these quantities can be very much 
influenced by strong electron correlation and metallic behavior, which are not exclusive properties 
of ferroelectric materials. In this paper we compare the Born effective charges of some prototypical 
ferroelectrics with those of magnetic, non-ferroelectric compounds using a novel, self-interaction 
free methodology that improves on the local-density approximation description of the electronic 
properties. We show that the inclusion of strong-correlation effects systematically reduces the 
size of the Born effective charges and the electron localization lengths. Furthermore we give an 
interpretation of the Born effective charges in terms of band energy structure and orbital occupations 
which can be used as a guideline to rationalize their values in the general case. 

PACS numbers: Valid Pacs 
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I. INTRODUCTION 

Following the implementation of the linear-response 
theory«i and the Berry phase approach^ to the electron 
polarization within the local-spin density approximation 
(LSD A), there has been tremendous progress in under- 
standing the properties of piezoelectric and ferroelectric 
materials from a microscopic standpoint. Indeed, spon- 
taneous polarization, piezoelectric constants, and Born 
effective charges (BEC), as well as full phonon spectra 
and potential energy surfaces can be calculated from first- 
principles giving thus a detailed picture of the electronic 
and structural behavior. 

The Born effective charges (BEC), i.e. the change in 
electron polarization upon ionic displacements, are the 
main subject of this paper. They play an important role 
in the physics of ferroelectric and even piezolectric ma- 
terials, since they contribute to the LO-TO phonon fre- 
quency splitting at T (which is characteristically large in 
ferroelectric materials^) and to the piezoelectric tensor^ 
Furthermore, the behavior of the BEC is highly non- 
trivial since it is related to the microscopic electron cur- 
rent produced in a system by a change of the atomic 
position. 4 This current may give a contribution to the 
BEC which is unrelated and additional to the static elec- 
tron charge carried by the ion during the displacement, 
so that the BEC can be much bigger than their static 
counterpart. In this case the BEC are usually referred to 
as anomalous. 

Interestingly, the BEC are found to be highly anoma- 
lous in the high-temperature (i.e. centrosymmetric) 
phase of many ferroelectric compounds. In particu- 
lar, they are anomalous for the whole family of ferro- 
electric perovskites^*2*2i2iiS and for the class of IV- VI 
chalcogenidesjii. Since both the BEC anomaly and the 
mechanism of the ferroelectric phase transition in most 
of the ferroelectric perovskites can be related to the hy- 
bridization of anion p and cation d orbitalsi 1 ^* 1 ^, the in- 



teresting hypothesis arises that a large value of the BEC 
might be an indicator of the ferroelectric materials, or 
materials which are on the verge of a ferroelectric transi- 
tion. However, according to recent calculations the BEC 
do not appear to posses such a predictive capability in 
general. On one hand they are found to be anomalous 
in a magnetic perovskite (CaMnOa 1 ^) which is far from 
ferroelectric. On the other hand they are non-anomalous 
in the ferroelectromagnet YMnOs, which belongs to the 
family of the hexagonal manganese oxides 1 ^* 1 ^* 1 ^, whose 
ferroelectric transition is driven not by chemical activity 
(e.g. anion-cation hybridization changes) but by struc- 
tural collapse and oxygen rotations^ 

So we see that the meaning of the BEC and origin of 
the anomalous part are far from being completely clari- 
fied. While the role of covalent bonding in favoring spon- 
taneous polarization and large BEC in ferroelectric per- 
ovskites is well understoodjZi2ii£ other aspects which can 
strongly influence the values of the BEC, like the metallic 
character and the strong electron correlation correlation, 
have been scarcely investigated so far. 

In an insulator on the verge of an insulating-to-metal 
transition, the BEC can diverge to anomalously large 
values, thus indicating an electron current which can 
freely flow through the system in response to the electric 
field associated with an ionic displacement. 1 - Even more 
interesting is the effect of strong electron-correlation. 
In a series of papers Ishihara and co-workers^ showed 
that in a two-band, Hubbard model, the introduction of 
strongly-correlated effects (namely the on-site Coulomb 
energy) induces a strong enhancement of the electron- 
phonon coupling which becomes divergent at the tran- 
sition boundary from the weakly-correlated (small U) 
to the strongly-correlated regime (large U). Resta and 
Sorella 20 showed that, in the same model, in the region 
of the transition boundary, the BEC are divergent and 
undergo a sign change (i.e. the electron current changes 
its direction flowing from the cation to the anion) for 
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large U. 

Later the same authors introduced the concept of elec- 
tron localization length^ (A) as a measure to discrimi- 
nate between the metallic regime (where A=oo) and the 
insulating regime. Within the same Hubbard model, they 
found that the BEC and A diverge simultaneously for the 
same value of U (i.e. the BEC diverge in the metallic 
state), 

Despite these predictions from model Hamiltonians, to 
our knowledge very few attempts have been made to cal- 
culate these properties for strongly-correlated materials 
from first-principles. The main reason is likely the ques- 
tionable reliabilty of LSDA calculations for describing 
the electronic properties of strongly localized, atomic- 
like, electron states. It is common knowledge, indeed, 
that the LSDA is often inaccurate in describing these 
systems However, there is increasing evidence that 
many of the LSDA failures are not due to the lack of gen- 
uine many-body electron correlation, but rather to the 
presence of the spurious self interaction (SI) within the 
LSDA functional^ which becomes particularly strong 
for spatially localized electron charges. The effects of the 
SI were first studied by Perdew and Zungen 2 ^ on atoms. 
Subsequently Svane and co-workers 24 implemented a self- 
interaction-corrected (SIC) scheme for extended systems, 
obtaining large improvements over the LSDA results for 
the electronic properties of transition metal oxides and 
high-T c superconductors. However their method was not 
simple and general enough to be extensively applied to 
a vast range of systems, thus the capabilities and the 
performance of the SIC remained largely untested. 

Recently we developed a self-interaction free pseudopo- 
tential scheme (pseudo-SIC) 2 ^ based on the idea (origi- 
nally conceived by Vogel and co-workers 2 ^) of express- 
ing the SI as a non-local, pseudopotential-like projec- 
tor. This approach has the advantage of being appli- 
cable to a vast generality of systems including semicon- 
ductors and metals, magnetic and non magnetic mate- 
rials, without requiring large modifications or a large 
increase of computing effort with respect to the usual 
LSDA. The pseudo-SIC is close in spirit to the popu- 
lar LDA+U2L2& since removing the SI is analogous to 
effectively mapping the on-site energy U into the one- 
electron picture, so that both weakly-correlated materials 
and Mott-insulators can be treated on the same foot. In 
Refl25l we showed that the method is accurate in describ- 
ing the band energies of III-V and II- VI semiconductors 
(e.g. GaN and ZnO), transition metal monoxides (NiO, 
MnO) and hexagonal ferroelectrics (YMJ1O3). 

In this paper we apply the pseudo-SIC to the study 
of the BEC for a range of different materials includ- 
ing a prototypical ferroelectric perovskite (BaTiOa), a 
ferroelectric with extremely high BEC (GeTe), a mag- 
netic perovskites (CaMnOa) and two magnetic monox- 
ides (MnO and NiO). This investigation aims to furnish 



a more complete and accurate picture of the BEC, and 
to spread light on the meaning of their anomalous values. 

We show that the BEC anomaly is not a property 
of ferroelectric materials, but is instead present in ionic 
insulators whenever the energy difference between filled 
anion and empty cation bands is particularly small (less 
than 3 eV), or in materials with nearly metallic behavior. 
Furthermore, we give evidence that the effect of strong 
correlation systematically reduces the value of the BEC 
and of the electron localization length. This decrease 
is particularly evident for the BEC which are descibed 
as highly anomalous within LDA (or LSDA). Thus, part 
of (but not all) the BEC anomaly found within LDA 
is a consequence of the unphysical self-interaction which 
overemphasizes the orbital hybridization and the sensi- 
tivity of the electron-phonon coupling to the atomic dis- 
placement. 

Finally, we propose an interpretation of the BEC in 
terms of occupation and capacity of the filled and empty 
orbitals, respectively, which allows a qualitative under- 
standing of the anomaly in a general case. 

The remainder of the paper is organized as follows: in 
Section |H] we briefly illustrate the pseudo-SIC scheme, 
and the formulation used for the BEC calculation. Sec- 
tion lllll is devoted to the presentation and the discussion 
of our results, and in Section[3we give our final remarks. 



II. FORMULATION 

A. Self-interaction corrected approach 

Here we briefly summarize the main features of the 
pseudo-SIC approach. A detailed discussion of the for- 
malism and the assumptions which are at the basis of 
this method can be found in Ref. 25. The formalism uses 
ultrasoft pseudopotentials (USPP™) which give optimal 
transferability and require low cut-off energies even for 
transition metal atoms. The Si-free, Kohn-Sham equa- 
tions are: 
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where Vpp is the non-local pseudopotential, the 

screening potential, and the SIC term. This is ex- 
pressed as a fully-non local, Kleinman-Bylander-type, 
pseudopotential projector: 
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Here i is a composite index which runs over angular 
quantum numbers and atomic positions, </>i(r) are atomic 
pseudowavefunctions, ^j^- c ( r ) are the SI potentials of 
the fully occupied spin-polarized atomic pseudocharges 
4>i(v) 2 , and pf are the corresponding occupation num- 
bers. The second term in Eq. [3 is the augmented con- 
tribution required if using USPP (r) and Q vlll (r) are 
projector functions and augmented charge desities^ re- 
spectively), whereas for ordinary norm-conserving pseu- 
dopotentials only the first term in Eq. [21 has to be con- 
sidered. 

All the ingredients needed in Eqs[21are atomic quanti- 
ties (thus they are calculated during the initialization), 
except the occupation numbers, pf , which must be evalu- 
ated self-consistently by projecting the Bloch state man- 
ifold onto the pscudoatomic basis: 



pr = E^« k i^){^iCk> 



(3) 



Again, the second term of Eq. [3] is only present 
for USPP. (The augmented charges q vil are integrals of 



B. Electron polarization and localization 

For the calculation of the BEC we employ the com- 
pact formulation developed by Resta 3 — , which is based 
on the concept of position operator. This has the advan- 
tage of introducing electronic polarization and electron 
localization length as two quantities steaming from the 
same source. Here we describe the formulation for spin- 
polarized systems and USPP which, to our knowledge, 
has not been presented previously. 

The source of the theory is the expectation value z a of 
the many-body phase operator: 



= (*o|e 



Ak„-R 



(4) 



Here a is the cartesian coordinate, R = YliLi r « the 
many-body position operator, \l/o the ground-state wave- 
function, and Ak Q a reciprocal-space vector of magni- 
tude: 
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Here L a and M are length and number of unit cells 
(that is number of grid points in reciprocal space) in the 
a direction, respectively. 

Resta showed that the expected value of the position 
operator, within periodic boundary conditions (PBC), 
must be calculated as: 



(r a ) 



NG, 
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where G a =27r/L Q , N is the number of electrons, and 
Im ln2: a is the phase of z a . From Eq. [5] the remarkable 
result follows that, within PBC, the electronic polariza- 
tion P is: 
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On the same footing, the electron localization, A Q , (i.e. 
the spread of the single-particle position operator) can be 



calculated from the absolute value of as: 



A* = (rl) - (r a ) 
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(8) 



Although Eqs[7| and [S] are valid for any system of in- 
teracting electrons, z 1 ^ is in practice calculated within 
the usual independent-particle approximation. Once the 
many-body wavefunction, is expressed as a sin- 

gle Slater determinant of Bloch wavcfunctions we have 
that^i 



z% = det Si 



k,mk' ' det S nk 



where, 
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These are the elements of a N T x NT = (]\^M 3 )x (n£m 3 ) 

matrix (N^ and are the number of spin-up electrons 
and occupied bands, respectively). It is easy to see that 
they are only non- vanishing for k' = k + Ak„ . Thus the 
determinant in Eq|5|can be decomposed into a product of 
M 3 determinants (i.e. one for each point of the reciprocal 
space grid) of an xN^-dimensional matrix: 



detSl k mk , = J] detMT >m (k,k+ Ak Q ), (11) 

k 

where, 
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MT im (k,k + Ak Q ) = (^ k |e- lAk - r |^ T „, k+Akci ).(12) 

A complication is caused by the use of USPP, since 
with the USPP method the Bloch states are not orthonor- 
mal, thus the expected value of the phase operator must 
be augmented by an additional term which restores the 
contribution of the augmented atomic charges, qy(r)'s. 
Thus, the following contribution M„~J should be added 

to the matrix M„ m given in EalT^I 

M^(k,k + Ak a ) = (<kl$) (^l< k+ A k „) x 

e -iAk„.R, J drqij ( T ) e -i^« r, (13) 

where R; are atomic positions. 

The calculation of Eqs. I11II12I and 1 131 is not demand- 
ing since the product over k-point space requires a dense 
grid only in the a-direction, whereas few k-points in the 
subspace orthogonal to a are generally sufficient to ob- 
tain a well-converged result for electron polarization and 
localization length along direction a. 

III. RESULTS 
A. Technicalities 

Our local-spin density calculations-- are performed 
within a plane-wave basis set and employ ultra- 
soft pseudopotentialsS For the LDA calculations, Ti 
(3s 2 ),(3p 6 ) semi-core electrons are included as valence 
charge in order to improve the transferability of their 
respective pseudopotentials. This is not necessary in 
pseudo-SIC since these states are much lower in energy 
(see Refill for a discussion about this point) . In all cases 
we employ two beta functions for each angular chan- 
nel. Also, non-linear core corrections have been used for 
all the atoms, and scalar-relativistic effects are included 
within the USPP for Ca, Ba, Ti, and Ge. 

We use cut-off energies equal to 30 Ryd for CaMn03 
and 40 Ryd for all the other compounds. For the calcu- 
lation of the Berry phase, 6x6 special k-point grids are 
used for the integration in the plane orthogonal to the 
polarization direction, and up to 40 k-points to integrate 
along the parallel direction. 

B. BaTi0 3 and GeTe 

In this Section we investigate two ferroelectric materi- 
als. BaTiC>3 is one of the simplest and best known ferro- 



electric perovskites, and the LDA successfully describes 
the energetics of the three ferroelectric phase transi- 
tions (tetragonal, orthorhombic and rhombohedral), pro- 
vided that the experimental volumes are used for the 
calculations. 32 The BEC values are also well-established 
by a long series of LDA calculations, and their large 
anomaly is assumed to be a major indication of ferro- 
electric instability. Also, the LDA calculations revealed 
that, despite the mainly ionic character of BaTiC>3, some 
p-d hybridization is present^ and the sensitivity of this 
hybridization to off-center atomic displacements is a fun- 
damental force driving the system to the ferroelectric 
stateiH Despite these successes, the comparison be- 
tween the LDA band-structure energies and the pho- 
toemission spectra^ 3 , show typical discrepancies that are 
worthy of investigation. 

In Figure IIII Bl (left panel) our calculated LDA band 
energies of cubic BaTiC>3 at the LDA volume are shown. 
The bands are well separated into groups, each with a 
single dominant orbital character. For the direct energy 
gap at r we obtain the value 1.94 eV (the absolute en- 
ergy gap is indirect between R and V points and equal 
to 1.84 eV), in agreement with other LDA calculations, 34 
while, as usual, the experiments indicate a much larger 
gap (3.2 e\4?£). Also, the spread of the O p manifold 
is ~ 4.5 eV, in contrast to an experimental value of 5.5 
eVi^ 3 . Finally, the binding energies of the more localized, 
lowest-lying states (9.9 eV for the Ba 5p bands, 16.7 eV 
for the O 2s bands) are underestimated with respect to 
the corresponding photoemission energies (12.2 eV and 
18.8 eV) 

We see in the right panel of Figure IIII Bl that the 
pseudo-SIC repairs the LDA deficiencies to a large ex- 
tent: The energy gap at T is 3.13 eV, the bandwidth of 
the O p manifold is ~ 6 eV, and the energies of Ba 5p 
and O 2s bands ~ 11.57 eV and 18 eV, respectively. 

There is a long-standing debated about the interpre- 
tation of the mismatch between the LDA (or LSDA) and 
the observed energy gap. Briefly, if the mismatch is at- 
tributed to the LDA inaccuracy, this may, in principle, 
affect any ground-state property. If otherwise the mis- 
match steams from a discontinuity of the exact energy 
functional (i.e. it depends on the DFT itself) it has no 
consequences on the ground-state properties. In fact, 
both these aspects contribute to the mismatch. In the 
present case, the agreement of pseudo-SIC and experi- 
mental energy gap is an indication that the LDA error is 
the dominant contribution. 

In Tabic [I] the calculated values for the BEC are 
shown. The most important result is the remarkable de- 
crease of Z* for both Ti and Ot, when calculated within 
pseudo-SIC, with respect to their LDA counterparts. The 
pseudo-SIC values can still be described as anomalous if 
compared with the formal ionic charges, but the anomaly 
(1.88) is reduced by ~ 40% with respect to the LDA. 
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FIG. 1: Band structure of cubic BaTiOs calculated within 
LDA and pseudo-SIC. 



TABLE I: Born effective charges, Z*, and localization lengths, 
A, (in Bohr) calculated within LDA and pseudo-SIC for cubic 
BaTiOs. Or and Op label oxygens on-top of Ti and planar 
with Ti, respectively. 





Z Ti 


z b T 


Op 


A 


formal 2 


4 


-2 


-2 




LDA 2.70 


7.01 


-5.54 


-2.05 


0.93 


pseudo-SIC 2.61 


5.88 


-4.43 


-2.03 


0.97 



This result may be expected since typically the SIC in- 
creases the removal energies and the space localization of 
the electron charges. This, in turn, reduces the overlap 
and the hybridization which is generally overemphasized 
within LDA. 

In Tabled the calculated localization lengths, A, are 
also reported. In a many-body framework A is an exact 
measure of the electron localization. In the one-electron 
approximation, it measures the average spatial localiza- 
tion of a single electron, thus it is a useful indicator of 
the overall metallic character of a system. The small 
value obtained for BaTiOs (less than 1 Bohr) is typical 
of highly ionic compounds (compare this, for example, 
with the values obtained for the series of III-V semicon- 
ductors in Ref©. Also, notice that the LDA value is 
smaller than the pseudo-SIC result. This is caused by the 
presence of Ti 3s, 3p semicore states in the LDA valence 
charge (these states are in the core for the SIC calcula- 



TABLE II: Orbital occupation numbers for cubic BaTiOs cal- 
culated within LDA and pseudo-SIC. The numbers show the 
total orbital character for individual atoms. Ba s and p orbital 
occupations are relative to the low-lying semi-core states. 



LDA 



pseudo-SIC 





Ba Ti 


O 


Ba 


Ti 


O 


s 


1.64 0.36 


1.60 


1.65 


0.38 


1.61 


p 


5.63 1.30 


4.82 


5.65 


1.26 


4.99 


d t 2s 


0.81 1.10 




0.72 


0.84 




d e 9 


0.42 1.16 




0.35 


1.08 





tion). Wc will see later that, for the same valence states, 
the localization length is always smaller in pseudo-SIC. 

The effect of the SIC on the orbital hybridizations can 
be appreciated in Table ILTl where the orbital occupation 
numbers are reported. In pseudo-SIC the material is 
more ionic and the hybridizations between Ti d and O 
p orbitals are lowered. Indeed, in pseudo-SIC the Ti e g , 
and t2 g , and the O p orbital occupations are closer to 
their respective nominal values. 

Now we present our results for GeTe t38i38 i This is an in- 
teresting ferroelectric since it possesses extremely anoma- 
lous BEG— and shows electronic properties markedly dif- 
ferent from those of the most common ferroelectric pre- 
rovskites. At high-temperature it has the rocksalt struc- 
ture, and below T c ~ 650 K 38 it undergoes a transition to 
a ferroelectric rhombohedral structure, with polarization 
parallel to the [111] axis. 

In Figure ITTTEfl the LSDA and pseudo-SIC band struc- 
tures of rocksalt GeTe are shown. There are no d states 
or other highly localized states close to the energy gap re- 
gion, and in fact the band structure shows fairly dispersed 
states both in the valence and in the conduction regions. 
Also, in this case LDA and pseudo-SIC give very similar- 
results. This is consequence of the large hybridizations 
of the bands around the energy gap, with s and p or- 
bital characters heavily mixed. When the bands are not 
dominated by a well-defined orbital character, or are not 
localized on a single atomic site, the SIC corrections over 
the LDA eigenvalues are small and mainly consist of a 
global shift in energy of the whole band manifold A 

Due to the large band dispersion, GeTe is a quasi- 
metal, with a direct energy gap at L=[l/2,l/2,l/2] 27r/ao 
equal to 0.43 eV in LDA and 0.20 eV in pseudo-SIC. The 
metallicity is consistent with the huge LDA BEC value 
previously reportediii Our results give Zrp e = -11.55 elec- 
trons within LDA and Zrp„ = -8.27 electrons in pseudo- 
SIC. Thus almost 35% of the LDA anomalous charge is 
actually due the SI error. The metallic character and the 
effect of the SIC can be quantitatively estimated from the 
values of the localization length. In LDA we obtain A = 
2.81 Bohr, a huge value in comparison to that the highly 
ionic BaTiOs. This indicates the large (nearly metal- 
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FIG. 2: Band structure of rock-salt GeTe calculated within 
LSDA and pseudo-SIC. 



lie) spatial extention of the valence electron charges. In 
pseudo-SIC we have A = 2.47 Bohr, (that is 23% smaller 
than the LDA value). The simultaneous decrease of A and 
BEC obtained within pseudo-SIC is consistent with the 
idea that in the metallic limit, A=oo, the BEC calculated 
within the centrosymmetric structure should diverge. 

To conclude this section we illustrate the mechanism of 
the anomalous currents in BaTi03 and GeTe by means 
of the simple scheme shown in Figure Consider first 
the upper panel which refers to BaTi03. For simplic- 
ity, only Ti and 0<r, which are the ions with the most 
anomalous BEC, are considered. Upon the Ti displace- 
ment (indicated by the arrow) a +4 ionic charge is trans- 
ported upward. Simultaneously, a further electron cur- 
rent flows from the above oxygen towards the Ti, so that 
the positive charge effectively carried by the Ti along the 
displacement is larger than +4. Of course, the amount 
of the anomalous contribution depends on the detail of 
the energy spectrum, but we can make some considera- 
tions based on the capacity of the orbitals available to 
exchange charge. If we assume that the current can only 
flow from the highest occupied O orbitals to the lowest 
unoccupied Ti orbitals^i, we can conclude by a simple 
criterion of orbital filling that the maximum value al- 
lowed for the anomalous contribution must be +8, that 
is (2s 2 ) (2p 2 ) oxygen electrons transferred to the 3d Ti 
orbitals, and than the largest BEC allowed according to 
this argument is 4+8 = +12 for Ti, and -2-8 = -10 for 
Ot- This criterion must hold for all the octahedrally- 
coordinated oxide perovskites with a d° B-site transition 
metal, and indeed first-principle calculations support this 
prediction. See for example Ref0 where a list of BEC for 
many ferroelectric perovskites is reported. 

We can test the validity of this model for the IV- VI 
chalcogenides whose reported BEC are very anomalousjii 
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FIG. 3: Schematic diagram of anomalous charge flux in 
BaTiOs and GeTe. Small arrows indicate atomic displace- 
ments, long arrows show the directions of electron currents. 
For BaTiOs the anomalous charge can only flow from filled 
O 2s, 2p to empty Ti 3d orbitals and amounts to a maximum 
of 8 electrons. For GeTe two channels are available, giving a 
total maximum current of 10 electrons. The calculated values 
are indeed below these maximum limits. 



the biggest of whom is indeed found in GcTc. In order 
to justify the large value of the anomalous part (9.55 in 
our LDA calculation) we should consider that there are 
two available channels of electron current (shown in the 
lower panel of Fig|3|): one involves the current from the 
filled Ge (4s) 2 to the empty Te (6s) 2 , (6p) 6 orbitals, and 
the other from the filled Te (5s) 2 , (5p) 6 to the empty 
Ge (4p) 6 , (5s) 2 orbitals. The sum of these two contribu- 
tions give 10 electrons as a maximum limit of anomalous 
charge, which is indeed consistent with the calculated 
value. 



CaMnO; 
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In this Section we report our results for the G-type^ 
antiferromagnetic CaMnOsj 4 ^ Among the manganese 
perovskites, CaMnOs is the only compound stable in 
cubic symmetry since the Mn 4+ ion has the three ma- 
jority d t 2g orbitals fully occupied, and therefore does 
not undergoes Jahn- Teller distortions. Due to its formal 
ionic configuration Ca 2+ Mn 4+ 2 + this is the magnetic, 
non-ferroelectric perovskite which is better suited to be 
compared with BaTiOs . 
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In a previous paper— devoted to the analysis of orbital 
hybridizations we reported the fact that, within LSDA, 
CaMnC>3 shows anomalous BEC, with values similar to 
that of BaTiC>3, and we argued that the BEC anomaly 
could be, in fact, a common characteristic of the whole 
family of perovskite oxides. However, the two materials 
show substantial differences in the electronic properties. 
BaTi03 is definitely more ionic in character, whereas 
CaMnOa show a larger band dispersion and hybridiza- 
tion. Furthermore, the small band gap obtained within 
LSDA for CaMn0 3 suggests that the large BEC could 
be due to an incipient insulator-metal transition. These 
facts motivated us to revisit the properties of CaMn03 
within pseudo-SIC. 

In Figure 21 we report the band structure calculations 
within LSDA and pseudo-SIC. Within pseudo-SIC the 
bands are flatter and spread over a wider energy range. 
In pseudo-SIC the valence band manifold of O 2p and 
Mn d t2g states spans a region of 8 eV below the valence 
band top (VBT), where the same band manifold is 7 eV 
wide in LSDA. In both pseudo-SIC and LSDA the O 2s 
state are placed ~ 18 eV below the VBT. Within LSDA 
the fundamental energy gap is only 0.42 eV and occurs 
between points A=[l/2, 1/2,0] 2?r/a (the VBT) and T, 
whereas within pseduo-SIC the gap is direct at L and 
equal to 1.88 eV. In both LSDA and pseudo-SIC the low- 
est group of conduction bands come from the majority e g 
and the minority t2 g orbitals. (The tz g bands are flatter 
and less hybridized than the e g bands since they do not 
point towards the oxigens.) 

Comparing these results with an x-ray absorption spec- 
troscopy study^i shows that the pseudo-SIC is in overall 
better agreement with experiment than LSDA. In par- 
ticular, the pseudo-SIC reproduces well the observed dis- 
tance between the average positions of the lowest conduc- 
tion and highest valence bands (~ 3 eV), which is instead 
strongly underestimated in LSDA (~ 1.5 eV). Even the 
observed width of the O 2p-Mn t2 9 occupied manifold (~ 
9 eV) is in better agreement with the pseudo-SIC value, 
whereas both LSDA and pseudo-SIC correctly describe 
the material as a charge-transfer insulator (i.e. U > A 
in the language of Hubbard model 44 ). 

On the basis of the energy band results, we expect 
that the BEC in pseudo-SIC should be markedly differ- 
ent from the LSDA values, and the results in Table IIIII 
confirm this expectation: the values for Zjyj n and Zq^ 

are reduced by ~ 25% in pseudo-SIC, altough they still 
remain substantially anomalous. It is thus evident that 
part, but not all, of the BEC anomaly calculated in LSDA 
is due to the overestimation of the p-d hybridization. 
Furthermore, it is clear that anomalously large values 
of the BEC can be, as in this case, absolutely unrelated 
to any notion of ferroelectric behavior. (It is helpful to 
remember, indeed, that the energetics of CaMnC>3, dis- 
favours any ferroelectric distortion—) 
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FIG. 4: Band energies for G-type AFM CaMn0 3 calculated 
within LSDA and pseudo-SIC. 



TABLE III: Born effective charges for CaMnOs calculated 
within LSDA and pseudo-SIC. For Mn we also report the 
individual spin-up {Z^^) an( I spin-down {Zf^ n ) electronic 
contributions. A is the localization length (in Bohr). 





%n 


z Mn 


z Mn 




Up 


A 


formal 2 


4 


-3 





-2 


-2 




LSDA 2.74 


7.60 


-2.57 


3.16 


-6.75 


-1.80 


1.42 


pseudo-SIC 2.87 


6.03 


-2.73 


1.76 


-5.36 


-1.77 


1.29 



Finally, notice that, as in BaTi03, the BEC for the 
A-site cation and Op are slightly anomalous. 

In Table IIIII we also report the localization length. 
Within pseudo-SIC A is smaller by ~ 20% than in LSDA, 
consistent with the reduction obtained for the BEC. 

Finally, a fundamental question arises: how does the 
spin-polarization affect the values of the BEC? The total 
(i.e. spin-integrated) BEC does not contain the necessary 
information to answer this question. Thus, in Table IIIII 
we show the electronic component of the Mn BEC de- 
composed into its spin-up and spin-down contributions 
(see Equation |5J), so that: 



I i 



z Mn - z Mn + Z Mn + Z Mn 



(14) 



with = 7. The result can be rationalized as fol- 



8 



lows: the majority contribution -Zj^ n is substantially 
non- anomalous and close to its formal value -3. This 
means that the filled t 2g shell moves almost rigidly with 

the Mn 7+ core. Instead, ^^ n ^ is highly anomalous. The 

tjL states, which are at the bottom of the conduction 
band, can accept at most 3 electrons of anomalous cur- 
rent flowing from the O onto the Mn. In LSDA the en- 
ergy separation between the center of the tjj fl band man- 
ifold and the top of the O p valence bands is very small 
(~ 1 eV), thus the calculated anomalous, spin-down cur- 
rent (3.16) matches this acceptance limit (it is actually 
actually higher due to other secondary flux channels). 
Within pseudo-SIC, instead, the energy separation is ~ 
3 eV, and the anomaly becomes drastically reduced to 
1.76 electrons, that is almost half the LSDA value. On 
the other hand, the BEC contribution from the majority 
channel in pseudo-SIC is almost the same as in LSDA. 

Notice that we do not show the separated spin-up and 
spin-down contributions for O since they are equal by 
symmetry. 

In conclusion, in a magnetic material a correct inter- 
pretation of the BEC values can be only achieved if spin- 
up and spin-down components are disentangled, since 
their contribution to the integrated value can be highly 
anisotropic. 



D. Transition- metal Monoxides 

In this section we consider two magnetic, non- 
ferroelectric materials which are typical examples of sys- 
tems for which the LSDA fails to give an accurate descrip- 
tion, i.e. the AFM [111] A-type MnO and NiQ. 45 i 46 i 47 
The band structures of these compounds within both 
LSDA and pseudo-SIC have been already presented in 
Refill thus they will not be shown here, instead, here 
we briefly summarize the main conclusions of that inves- 
tigation: Within LSDA the fundamental energy gaps (0.9 
eV for MnO and 0.4 eV for Ni) are much smaller than the 
experimental values, and are of Mott-Hubbard type (i.e. 
the lowest energy excitation occurs between filled and 
empty d bands). Furthermore, the local magnetic mo- 
ments (4.42 [is and 1.11 [is for Mn and Ni, respectively) 
underestimate the measured values. The pseudo-SIC re- 
stores a good overall agreement between band energies 
and photoemission spectra: the energy gaps are 3.98 eV 
for MnO and 3.89 for NiO, the magnetic moments 4.71 
fis for Mn and 1.77 /is for Ni, and the character of the 
energy gap is in the intermediate charge-transfer Mott- 
Hubbard regime, in agreement with experiments. 

In Table llVl the calculated values for the BEC and A 
within LSDA and pseudo-SIC are shown. We first discuss 
MnO. Perhaps surprisingly after the strong differences 
described above for the band energy structures, the BEC 



TABLE IV: Born effective charges for MnO and NiO calcu- 
lated within LSDA and pseudo-SIC. The individual spin-up 
and spin-down electronic contributions are also shown. A is 
the localization length (in Bohr). 



MnO 


%n 


rvd T 

Z Mn 


iyel ■ 

Z Mn 


A 


formal 


2 


-5 







LSDA 


2.31 


-5.13 


0.44 


1.36 


pseudo-SIC 


2.29 


-4.96 


0.25 


1.19 


NiO 


z m 


ryel T 

Z Ni 


rye.1 I 

z Ni 


A 


formal 


2 


-5 


-3 




LSDA 


2.31 


-4.47 


-3.22 


1.18 


pseudo-SIC 


2.04 


-4.87 


-3.09 


0.76 



are found to be similar and substantially non anoma- 
lous within LSDA and pseudo-SIC (our LSDA result is in 
perfect agreement with an earlier linear-augmented plane 
waves calculation^) . This result can be understood from 
the band structures shown in Ref|25l MnO has, nomi- 
nally, five fully occupied majority d bands and five empty 
minority d bands. In LSDA these two manifolds are sep- 
arated by Hund's rule, whereas in pseudo-SIC the on-site 
Coulomb energy is dominant and the d-d energy gap in- 
creases dramatically. However, for this material the pa- 
rameter which controls the anomalous charge is not the 
fundamental gap, but the energy separation between the 
highest occupied O p bands and the lowest unoccupied 
Mn d band. This is large (~ 4 eV) and similar in both 
LSDA and pseudo-SIC, thus the non-anomalous value for 
the BEC are justified. 

Now consider the BEC decomposition within ionic and 
electronic, up-spin and down-spin components (Equation 
I14f) . The up-spin and down-spin formal charges are -5 
and 0. We see that in LSDA both up-spin and down-spin 
electronic contributions are slightly anomalous, but the 
anomalies in the two channels have opposite signs (that is 
the currents flow in opposite directions) thus they cancel 
each other in the integrated value. The pseudo-SIC, in 
fact, reduces the anomalies in both channels, but, due to 
the current cancellations, this is not revealed in the total 
BEC which is indeed very close to the LSDA value. 

For NiO we can repeat most of the considerations ex- 
pressed for MnO. Now the energy gap in LSDA lies in the 
minority d band manifold and is due to the crystal field 
splitting which is only ~ 1 eV. This value underestimates 
by far the experimental value. However, as in MnO the 
p-d energy separation is ~ 4 eV and therefore no large 
anomalous charges contribute to the BEC. From the up- 
spin and down-spin decomposition (now Z^V- = 10) we 
see that there is also a cancellation between the corre- 
sponding anomalies. The cancellation is almost complete 
within pseudo-SIC, thus the BEC is almost equal to the 
formal value. 
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Finally, within pseudo-SIC we obtain a significant de- 
crease of localization lengths for both these compounds 
with respect to the LSDA. 



IV. DISCUSSION 

Finally, we summarize the main findings of our analy- 
sis. We have seen that including strong-correlation effects 
within the band energy structure of strongly-correlated 
materials sensitively affects the values of the BEC and the 
electron localization lengths. With respect to the LDA 
or LSDA results, the pseudo-SIC systematically reduces 
both BEC and A for all the investigated materials. This 
effect is particularly evident for materials with highly 
anomalous BEC (i.e. BaTi03 and CaMnOa), and can 
be understood as resulting from the fact that the LSDA 
overemphasizes the orbital hybridizations and therefore 
the electron hopping. The decrease of electron local- 
ization lengths is also expected since the SIC enhances 
the binding energies and therefore the spatial localiza- 
tion of the electron charges. However, we point out that 
even within pseudo-SIC the anomalous contribution to 
the BEC is still substantial. 

Our investigation reveals that the BEC are controlled 
by two fundamental parameters: one is the the energy 
separation between the highest occupied anionic bands 
and the lowest empty cationic bands. This criterion 
explains the trend of the BEC calculated for BaTiOs 
and CaMnOs (where the separation is less than 2 eV 
in LSDA and ~ 3 eV in pseudo-SIC) and the transition- 
metal mono-oxides, where it is ~ 4 eV in both LSDA and 
pseudo-SIC. The other important element is the metal- 
lic character of the system. It is very sound that in the 
metallic limit the BEC must be divergent, and the elec- 
tron localization length gives us the possibility of evalu- 
ating this limit quantitatively. GeTe is a perfect material 
for such an analysis, since it is a ferroelectric with highly 
anomalous BEC and nearly metallic behavior. We found 
the GeTe energy gap small in both LDA and pseudo-SIC, 
but in pseudo-SIC both the localization length and the 
BEC are markedly smaller. 

Finally, we should comment on the results of the two- 
band Hubbard model calculationsjii^Si since our results 
are in contrast with these works, which predict that the 
introduction of the on-site Coulomb energy, U, should 
increase the BEC with respect to a single-particle scheme, 
and the BEC should undergo a sign change corresponding 
to the charge-transfer to Mott- Hubbard transition. 

The disagreement is due to the fact that in the Hub- 
bard model it is assumed that the on-site Coulomb energy 
of the anion favors a ground state with mixed anion- 
cation occupancy, instead of a purely anionic ground 
state. This causes, at large U, the transition from the 
charge-transfer to the Mott-Hubbard regime, and also 



produces a sign change in the BEC since now the anoma- 
lous current can flow from the filled cation to the empty 
anion. 

In density fuctional calculations what happens is quite 
the opposite: with the on-site Coulomb interaction on the 
ionic site turned on, a larger energy cost is associated 
with the fluctuation of the respective occupancy. The 
tendency to fluctuation is reduced by lowering (not rais- 
ing) the energy of the filled anion bands. Thus, strong 
correlation effects further favor charge-transfer ground 
states over Mott-Hubbard regimes. This is in complete 
agreement with the photoemission spectroscopy data^^i 
for transition-metal and perovskite oxides, as well as 
with the LDA+U results .SLiMS In conclusion, we think 
that the two-band Hubbard model is not adequate to de- 
scribe materials like the ones studied in this paper, which 
present a complicated cohexistence of charge localization 
and hybridization. 



V. CONCLUSIONS 

In this paper we revisited the behavior of the BEC 
by means of our recently proposed pseudo-SIC method 
which set the LSDA density functional free of the spu- 
rious SI, and thus describes the electronic properties of 
weakly as well as strongly correlated materials efficiently. 
We investigated two ferroelectric materials with highly 
anomalous BEC, and some magnetic materials with both 
anomalous and non-anomalous BEC. 

Overall, we find that the anomalous part is system- 
atically reduced (but not suppressed) by the SIC. The 
electron localization lengths is also reduced by the SIC, 
consistently with the idea that in the metallic limit the 
BEC must be divergent. 

On the basis of these considerations, what can be ar- 
gued about the BEC's ability to predict a ferroelectric 
instability? According to our results, there is no evi- 
dence to claim such a role for the BEC, since their val- 
ues are related to features of the band energy spectrum 
which are not specifically properties of a ferroelectric ma- 
terial. We may argue that a compound with large BEC, 
if ferroeoectric, will have a large spontaneous polariza- 
tion, since this represents the integrated value of the BEC 
over the ferroelectric displacement^ but the BEC value 
in itself is not indicative of, or specific to a chemical en- 
vironment particularly favorable for non-ferroelectric-to- 
ferroelectric phase transition. 
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